\(\int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx\) [46]

Optimal result

Integrand size = 36, antiderivative size = 153 \[ \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx=-\frac {a^3 (A+6 B) x}{c^3}+\frac {a^3 (A+6 B) \cos (e+f x)}{c^3 f}+\frac {a^3 (A+B) c^3 \cos ^7(e+f x)}{5 f (c-c \sin (e+f x))^6}-\frac {2 a^3 (A+6 B) c \cos ^5(e+f x)}{15 f (c-c \sin (e+f x))^4}+\frac {2 a^3 (A+6 B) c^3 \cos ^3(e+f x)}{3 f \left (c^3-c^3 \sin (e+f x)\right )^2} \]

[Out]

-a^3*(A+6*B)*x/c^3+a^3*(A+6*B)*cos(f*x+e)/c^3/f+1/5*a^3*(A+B)*c^3*cos(f*x+e)^7/f/(c-c*sin(f*x+e))^6-2/15*a^3*(
A+6*B)*c*cos(f*x+e)^5/f/(c-c*sin(f*x+e))^4+2/3*a^3*(A+6*B)*c^3*cos(f*x+e)^3/f/(c^3-c^3*sin(f*x+e))^2

Rubi [A] (verified)

Time = 0.25 (sec) , antiderivative size = 153, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.139, Rules used = {3046, 2938, 2759, 2761, 8} \[ \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx=\frac {a^3 (A+6 B) \cos (e+f x)}{c^3 f}+\frac {a^3 c^3 (A+B) \cos ^7(e+f x)}{5 f (c-c \sin (e+f x))^6}+\frac {2 a^3 c^3 (A+6 B) \cos ^3(e+f x)}{3 f \left (c^3-c^3 \sin (e+f x)\right )^2}-\frac {a^3 x (A+6 B)}{c^3}-\frac {2 a^3 c (A+6 B) \cos ^5(e+f x)}{15 f (c-c \sin (e+f x))^4} \]

[In]

Int[((a + a*Sin[e + f*x])^3*(A + B*Sin[e + f*x]))/(c - c*Sin[e + f*x])^3,x]

[Out]

-((a^3*(A + 6*B)*x)/c^3) + (a^3*(A + 6*B)*Cos[e + f*x])/(c^3*f) + (a^3*(A + B)*c^3*Cos[e + f*x]^7)/(5*f*(c - c
*Sin[e + f*x])^6) - (2*a^3*(A + 6*B)*c*Cos[e + f*x]^5)/(15*f*(c - c*Sin[e + f*x])^4) + (2*a^3*(A + 6*B)*c^3*Co
s[e + f*x]^3)/(3*f*(c^3 - c^3*Sin[e + f*x])^2)

Rule 8

Int[a_, x_Symbol] :> Simp[a*x, x] /; FreeQ[a, x]

Rule 2759

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> Simp[2*g*(g
*Cos[e + f*x])^(p - 1)*((a + b*Sin[e + f*x])^(m + 1)/(b*f*(2*m + p + 1))), x] + Dist[g^2*((p - 1)/(b^2*(2*m +
p + 1))), Int[(g*Cos[e + f*x])^(p - 2)*(a + b*Sin[e + f*x])^(m + 2), x], x] /; FreeQ[{a, b, e, f, g}, x] && Eq
Q[a^2 - b^2, 0] && LeQ[m, -2] && GtQ[p, 1] && NeQ[2*m + p + 1, 0] &&  !ILtQ[m + p + 1, 0] && IntegersQ[2*m, 2*
p]

Rule 2761

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)/((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]), x_Symbol] :> Simp[g*((g*Cos[e
+ f*x])^(p - 1)/(b*f*(p - 1))), x] + Dist[g^2/a, Int[(g*Cos[e + f*x])^(p - 2), x], x] /; FreeQ[{a, b, e, f, g
}, x] && EqQ[a^2 - b^2, 0] && GtQ[p, 1] && IntegerQ[2*p]

Rule 2938

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.)
+ (f_.)*(x_)]), x_Symbol] :> Simp[(b*c - a*d)*(g*Cos[e + f*x])^(p + 1)*((a + b*Sin[e + f*x])^m/(a*f*g*(2*m + p
+ 1))), x] + Dist[(a*d*m + b*c*(m + p + 1))/(a*b*(2*m + p + 1)), Int[(g*Cos[e + f*x])^p*(a + b*Sin[e + f*x])^
(m + 1), x], x] /; FreeQ[{a, b, c, d, e, f, g, m, p}, x] && EqQ[a^2 - b^2, 0] && (LtQ[m, -1] || ILtQ[Simplify[
m + p], 0]) && NeQ[2*m + p + 1, 0]

Rule 3046

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Dist[a^m*c^m, Int[Cos[e + f*x]^(2*m)*(c + d*Sin[e + f*x])^(n - m)*(A + B
*Sin[e + f*x]), x], x] /; FreeQ[{a, b, c, d, e, f, A, B, n}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && I
ntegerQ[m] &&  !(IntegerQ[n] && ((LtQ[m, 0] && GtQ[n, 0]) || LtQ[0, n, m] || LtQ[m, n, 0]))

Rubi steps \begin{align*} \text {integral}& = \left (a^3 c^3\right ) \int \frac {\cos ^6(e+f x) (A+B \sin (e+f x))}{(c-c \sin (e+f x))^6} \, dx \\ & = \frac {a^3 (A+B) c^3 \cos ^7(e+f x)}{5 f (c-c \sin (e+f x))^6}-\frac {1}{5} \left (a^3 (A+6 B) c^2\right ) \int \frac {\cos ^6(e+f x)}{(c-c \sin (e+f x))^5} \, dx \\ & = \frac {a^3 (A+B) c^3 \cos ^7(e+f x)}{5 f (c-c \sin (e+f x))^6}-\frac {2 a^3 (A+6 B) c \cos ^5(e+f x)}{15 f (c-c \sin (e+f x))^4}+\frac {1}{3} \left (a^3 (A+6 B)\right ) \int \frac {\cos ^4(e+f x)}{(c-c \sin (e+f x))^3} \, dx \\ & = \frac {a^3 (A+B) c^3 \cos ^7(e+f x)}{5 f (c-c \sin (e+f x))^6}-\frac {2 a^3 (A+6 B) c \cos ^5(e+f x)}{15 f (c-c \sin (e+f x))^4}+\frac {2 a^3 (A+6 B) \cos ^3(e+f x)}{3 c f (c-c \sin (e+f x))^2}-\frac {\left (a^3 (A+6 B)\right ) \int \frac {\cos ^2(e+f x)}{c-c \sin (e+f x)} \, dx}{c^2} \\ & = \frac {a^3 (A+6 B) \cos (e+f x)}{c^3 f}+\frac {a^3 (A+B) c^3 \cos ^7(e+f x)}{5 f (c-c \sin (e+f x))^6}-\frac {2 a^3 (A+6 B) c \cos ^5(e+f x)}{15 f (c-c \sin (e+f x))^4}+\frac {2 a^3 (A+6 B) \cos ^3(e+f x)}{3 c f (c-c \sin (e+f x))^2}-\frac {\left (a^3 (A+6 B)\right ) \int 1 \, dx}{c^3} \\ & = -\frac {a^3 (A+6 B) x}{c^3}+\frac {a^3 (A+6 B) \cos (e+f x)}{c^3 f}+\frac {a^3 (A+B) c^3 \cos ^7(e+f x)}{5 f (c-c \sin (e+f x))^6}-\frac {2 a^3 (A+6 B) c \cos ^5(e+f x)}{15 f (c-c \sin (e+f x))^4}+\frac {2 a^3 (A+6 B) \cos ^3(e+f x)}{3 c f (c-c \sin (e+f x))^2} \\ \end{align*}

Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(316\) vs. \(2(153)=306\).

Time = 11.41 (sec) , antiderivative size = 316, normalized size of antiderivative = 2.07 \[ \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx=\frac {a^3 \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right ) \left (24 (A+B) \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )-4 (11 A+21 B) \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^3-15 (A+6 B) (e+f x) \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^5+15 B \cos (e+f x) \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^5+48 (A+B) \sin \left (\frac {1}{2} (e+f x)\right )-8 (11 A+21 B) \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^2 \sin \left (\frac {1}{2} (e+f x)\right )+4 (23 A+93 B) \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right )^4 \sin \left (\frac {1}{2} (e+f x)\right )\right ) (1+\sin (e+f x))^3}{15 f \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^6 (c-c \sin (e+f x))^3} \]

[In]

Integrate[((a + a*Sin[e + f*x])^3*(A + B*Sin[e + f*x]))/(c - c*Sin[e + f*x])^3,x]

[Out]

(a^3*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])*(24*(A + B)*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2]) - 4*(11*A + 21*B)
*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])^3 - 15*(A + 6*B)*(e + f*x)*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])^5 + 15
*B*Cos[e + f*x]*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])^5 + 48*(A + B)*Sin[(e + f*x)/2] - 8*(11*A + 21*B)*(Cos[(
e + f*x)/2] - Sin[(e + f*x)/2])^2*Sin[(e + f*x)/2] + 4*(23*A + 93*B)*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])^4*S
in[(e + f*x)/2])*(1 + Sin[e + f*x])^3)/(15*f*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^6*(c - c*Sin[e + f*x])^3)

Maple [A] (verified)

Time = 0.92 (sec) , antiderivative size = 157, normalized size of antiderivative = 1.03

method	result	size
derivativedivides	\(\frac {2 a^{3} \left (\frac {B}{1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )}-\left (A +6 B \right ) \arctan \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )\right )-\frac {2 A +6 B}{\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1}-\frac {8 A -8 B}{2 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{2}}-\frac {32 A +32 B}{5 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{5}}-\frac {40 A +24 B}{3 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{3}}-\frac {64 A +64 B}{4 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{4}}\right )}{f \,c^{3}}\)	\(157\)
default	\(\frac {2 a^{3} \left (\frac {B}{1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )}-\left (A +6 B \right ) \arctan \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )\right )-\frac {2 A +6 B}{\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1}-\frac {8 A -8 B}{2 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{2}}-\frac {32 A +32 B}{5 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{5}}-\frac {40 A +24 B}{3 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{3}}-\frac {64 A +64 B}{4 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{4}}\right )}{f \,c^{3}}\)	\(157\)
risch	\(-\frac {a^{3} x A}{c^{3}}-\frac {6 a^{3} x B}{c^{3}}+\frac {B \,a^{3} {\mathrm e}^{i \left (f x +e \right )}}{2 c^{3} f}+\frac {B \,a^{3} {\mathrm e}^{-i \left (f x +e \right )}}{2 c^{3} f}+\frac {-\frac {112 A \,a^{3} {\mathrm e}^{2 i \left (f x +e \right )}}{3}-24 i A \,a^{3} {\mathrm e}^{3 i \left (f x +e \right )}+\frac {56 i A \,a^{3} {\mathrm e}^{i \left (f x +e \right )}}{3}+12 A \,a^{3} {\mathrm e}^{4 i \left (f x +e \right )}-144 B \,a^{3} {\mathrm e}^{2 i \left (f x +e \right )}-104 i B \,a^{3} {\mathrm e}^{3 i \left (f x +e \right )}+88 i B \,a^{3} {\mathrm e}^{i \left (f x +e \right )}+36 B \,a^{3} {\mathrm e}^{4 i \left (f x +e \right )}+\frac {92 A \,a^{3}}{15}+\frac {124 B \,a^{3}}{5}}{\left ({\mathrm e}^{i \left (f x +e \right )}-i\right )^{5} f \,c^{3}}\)	\(223\)
parallelrisch	\(-\frac {\left (\left (\left (\frac {233}{2}-60 f x \right ) B -10 f x A +24 A \right ) \cos \left (\frac {f x}{2}+\frac {e}{2}\right )+\left (\left (-\frac {33}{2}+30 f x \right ) B +5 f x A -\frac {16 A}{3}\right ) \cos \left (\frac {3 f x}{2}+\frac {3 e}{2}\right )+\left (f x A +6 f x B -\frac {243}{10} B -\frac {24}{5} A \right ) \cos \left (\frac {5 f x}{2}+\frac {5 e}{2}\right )+\left (\left (-\frac {143}{2}+60 f x \right ) B -\frac {32 A}{3}+10 f x A \right ) \sin \left (\frac {f x}{2}+\frac {e}{2}\right )+\left (\left (-\frac {155}{2}+30 f x \right ) B +5 f x A -12 A \right ) \sin \left (\frac {3 f x}{2}+\frac {3 e}{2}\right )+\left (\left (-6 f x -\frac {11}{2}\right ) B -f x A -\frac {4 A}{3}\right ) \sin \left (\frac {5 f x}{2}+\frac {5 e}{2}\right )-\frac {B \left (\cos \left (\frac {7 f x}{2}+\frac {7 e}{2}\right )-\sin \left (\frac {7 f x}{2}+\frac {7 e}{2}\right )\right )}{2}\right ) a^{3}}{f \,c^{3} \left (-10 \cos \left (\frac {f x}{2}+\frac {e}{2}\right )+5 \cos \left (\frac {3 f x}{2}+\frac {3 e}{2}\right )-\sin \left (\frac {5 f x}{2}+\frac {5 e}{2}\right )+5 \sin \left (\frac {3 f x}{2}+\frac {3 e}{2}\right )+10 \sin \left (\frac {f x}{2}+\frac {e}{2}\right )+\cos \left (\frac {5 f x}{2}+\frac {5 e}{2}\right )\right )}\)	\(265\)
norman	\(\frac {-\frac {5 a^{3} \left (A +6 B \right ) x \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{c}-\frac {30 a^{3} \left (A +6 B \right ) x \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}+\frac {5 a^{3} \left (A +6 B \right ) x \left (\tan ^{12}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}-\frac {14 a^{3} \left (A +6 B \right ) x \left (\tan ^{11}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}+\frac {a^{3} \left (A +6 B \right ) x}{c}-\frac {a^{3} \left (A +6 B \right ) x \left (\tan ^{13}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}+\frac {\left (40 A \,a^{3}+246 B \,a^{3}\right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{3 c f}+\frac {14 a^{3} \left (A +6 B \right ) x \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}-\frac {84 a^{3} \left (A +6 B \right ) x \left (\tan ^{7}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}+\frac {71 a^{3} \left (A +6 B \right ) x \left (\tan ^{8}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}+\frac {84 a^{3} \left (A +6 B \right ) x \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}+\frac {51 a^{3} \left (A +6 B \right ) x \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}-\frac {71 a^{3} \left (A +6 B \right ) x \left (\tan ^{5}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}+\frac {30 a^{3} \left (A +6 B \right ) x \left (\tan ^{10}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}-\frac {51 a^{3} \left (A +6 B \right ) x \left (\tan ^{9}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c}-\frac {52 A \,a^{3}+282 B \,a^{3}}{15 c f}-\frac {\left (2012 A \,a^{3}+8802 B \,a^{3}\right ) \left (\tan ^{8}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{15 c f}+\frac {\left (136 A \,a^{3}+966 B \,a^{3}\right ) \left (\tan ^{9}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 c f}-\frac {2 \left (64 A \,a^{3}+255 B \,a^{3}\right ) \left (\tan ^{10}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 c f}+\frac {2 \left (4 A \,a^{3}+29 B \,a^{3}\right ) \left (\tan ^{11}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c f}-\frac {\left (4 A \,a^{3}+12 B \,a^{3}\right ) \left (\tan ^{12}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c f}-\frac {2 \left (304 A \,a^{3}+1539 B \,a^{3}\right ) \left (\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{15 c f}+\frac {2 \left (92 A \,a^{3}+591 B \,a^{3}\right ) \left (\tan ^{3}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 c f}-\frac {\left (1972 A \,a^{3}+9552 B \,a^{3}\right ) \left (\tan ^{4}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{15 c f}+\frac {\left (112 A \,a^{3}+748 B \,a^{3}\right ) \left (\tan ^{5}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{c f}-\frac {4 \left (712 A \,a^{3}+3297 B \,a^{3}\right ) \left (\tan ^{6}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{15 c f}+\frac {4 \left (76 A \,a^{3}+525 B \,a^{3}\right ) \left (\tan ^{7}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )}{3 c f}}{\left (1+\tan ^{2}\left (\frac {f x}{2}+\frac {e}{2}\right )\right )^{4} c^{2} \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{5}}\)	\(772\)

[In]

int((a+a*sin(f*x+e))^3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^3,x,method=_RETURNVERBOSE)

[Out]

2/f*a^3/c^3*(B/(1+tan(1/2*f*x+1/2*e)^2)-(A+6*B)*arctan(tan(1/2*f*x+1/2*e))-(2*A+6*B)/(tan(1/2*f*x+1/2*e)-1)-1/
2*(8*A-8*B)/(tan(1/2*f*x+1/2*e)-1)^2-1/5*(32*A+32*B)/(tan(1/2*f*x+1/2*e)-1)^5-1/3*(40*A+24*B)/(tan(1/2*f*x+1/2
*e)-1)^3-1/4*(64*A+64*B)/(tan(1/2*f*x+1/2*e)-1)^4)

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 337 vs. \(2 (150) = 300\).

Time = 0.26 (sec) , antiderivative size = 337, normalized size of antiderivative = 2.20 \[ \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx=\frac {15 \, B a^{3} \cos \left (f x + e\right )^{4} + 60 \, {\left (A + 6 \, B\right )} a^{3} f x - 24 \, {\left (A + B\right )} a^{3} - {\left (15 \, {\left (A + 6 \, B\right )} a^{3} f x - {\left (46 \, A + 231 \, B\right )} a^{3}\right )} \cos \left (f x + e\right )^{3} - {\left (45 \, {\left (A + 6 \, B\right )} a^{3} f x + 2 \, {\left (A + 66 \, B\right )} a^{3}\right )} \cos \left (f x + e\right )^{2} + 6 \, {\left (5 \, {\left (A + 6 \, B\right )} a^{3} f x - 2 \, {\left (6 \, A + 31 \, B\right )} a^{3}\right )} \cos \left (f x + e\right ) - {\left (15 \, B a^{3} \cos \left (f x + e\right )^{3} + 60 \, {\left (A + 6 \, B\right )} a^{3} f x + 24 \, {\left (A + B\right )} a^{3} - {\left (15 \, {\left (A + 6 \, B\right )} a^{3} f x + 2 \, {\left (23 \, A + 108 \, B\right )} a^{3}\right )} \cos \left (f x + e\right )^{2} + 6 \, {\left (5 \, {\left (A + 6 \, B\right )} a^{3} f x - 2 \, {\left (4 \, A + 29 \, B\right )} a^{3}\right )} \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )}{15 \, {\left (c^{3} f \cos \left (f x + e\right )^{3} + 3 \, c^{3} f \cos \left (f x + e\right )^{2} - 2 \, c^{3} f \cos \left (f x + e\right ) - 4 \, c^{3} f - {\left (c^{3} f \cos \left (f x + e\right )^{2} - 2 \, c^{3} f \cos \left (f x + e\right ) - 4 \, c^{3} f\right )} \sin \left (f x + e\right )\right )}} \]

[In]

integrate((a+a*sin(f*x+e))^3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^3,x, algorithm="fricas")

[Out]

1/15*(15*B*a^3*cos(f*x + e)^4 + 60*(A + 6*B)*a^3*f*x - 24*(A + B)*a^3 - (15*(A + 6*B)*a^3*f*x - (46*A + 231*B)
*a^3)*cos(f*x + e)^3 - (45*(A + 6*B)*a^3*f*x + 2*(A + 66*B)*a^3)*cos(f*x + e)^2 + 6*(5*(A + 6*B)*a^3*f*x - 2*(
6*A + 31*B)*a^3)*cos(f*x + e) - (15*B*a^3*cos(f*x + e)^3 + 60*(A + 6*B)*a^3*f*x + 24*(A + B)*a^3 - (15*(A + 6*
B)*a^3*f*x + 2*(23*A + 108*B)*a^3)*cos(f*x + e)^2 + 6*(5*(A + 6*B)*a^3*f*x - 2*(4*A + 29*B)*a^3)*cos(f*x + e))
*sin(f*x + e))/(c^3*f*cos(f*x + e)^3 + 3*c^3*f*cos(f*x + e)^2 - 2*c^3*f*cos(f*x + e) - 4*c^3*f - (c^3*f*cos(f*
x + e)^2 - 2*c^3*f*cos(f*x + e) - 4*c^3*f)*sin(f*x + e))

Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 4665 vs. \(2 (143) = 286\).

Time = 13.60 (sec) , antiderivative size = 4665, normalized size of antiderivative = 30.49 \[ \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx=\text {Too large to display} \]

[In]

integrate((a+a*sin(f*x+e))**3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))**3,x)

[Out]

Piecewise((-15*A*a**3*f*x*tan(e/2 + f*x/2)**7/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 +
165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f
*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) + 75*A*a**3*f*x*tan(e/2 + f*x/2)**6/(15*c**3*f*
tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*
x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c*
*3*f) - 165*A*a**3*f*x*tan(e/2 + f*x/2)**5/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 16
5*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*ta
n(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) + 225*A*a**3*f*x*tan(e/2 + f*x/2)**4/(15*c**3*f*ta
n(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/
2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3
*f) - 225*A*a**3*f*x*tan(e/2 + f*x/2)**3/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*
c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(
e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) + 165*A*a**3*f*x*tan(e/2 + f*x/2)**2/(15*c**3*f*tan(
e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)
**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f
) - 75*A*a**3*f*x*tan(e/2 + f*x/2)/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f
*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 +
f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) + 15*A*a**3*f*x/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f
*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 +
f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 60*A*a**3*tan(e/2 + f*x
/2)**6/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c
**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/
2 + f*x/2) - 15*c**3*f) + 120*A*a**3*tan(e/2 + f*x/2)**5/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 +
f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 -
165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 460*A*a**3*tan(e/2 + f*x/2)**4/(15
*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(
e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2)
- 15*c**3*f) + 320*A*a**3*tan(e/2 + f*x/2)**3/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6
+ 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*
f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 452*A*a**3*tan(e/2 + f*x/2)**2/(15*c**3*f*ta
n(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/
2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3
*f) + 200*A*a**3*tan(e/2 + f*x/2)/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*
tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f
*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 52*A*a**3/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(
e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2
)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 90*B*a**3*f*x*tan(e/2 + f*x/
2)**7/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c*
*3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2
+ f*x/2) - 15*c**3*f) + 450*B*a**3*f*x*tan(e/2 + f*x/2)**6/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2
+ f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**
3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 990*B*a**3*f*x*tan(e/2 + f*x/2)
**5/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3
*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 +
f*x/2) - 15*c**3*f) + 1350*B*a**3*f*x*tan(e/2 + f*x/2)**4/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2
+ f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3
- 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 1350*B*a**3*f*x*tan(e/2 + f*x/2)
**3/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3
*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 +
f*x/2) - 15*c**3*f) + 990*B*a**3*f*x*tan(e/2 + f*x/2)**2/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 +
f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3
- 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 450*B*a**3*f*x*tan(e/2 + f*x/2)/(
15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*ta
n(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/
2) - 15*c**3*f) + 90*B*a**3*f*x/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*ta
n(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x
/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 180*B*a**3*tan(e/2 + f*x/2)**6/(15*c**3*f*tan(e/2 + f*x/2)*
*7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**
3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) + 870*B*a**
3*tan(e/2 + f*x/2)**5/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*
x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 7
5*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) - 2010*B*a**3*tan(e/2 + f*x/2)**4/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c
**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e
/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) + 2220*B*a**3*tan(e/
2 + f*x/2)**3/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5
- 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f
*tan(e/2 + f*x/2) - 15*c**3*f) - 2232*B*a**3*tan(e/2 + f*x/2)**2/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*ta
n(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x
/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f) + 1230*B*a**3*tan(e/2 + f*x/
2)/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*tan(e/2 + f*x/2)**5 - 225*c**3*
f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*x/2)**2 + 75*c**3*f*tan(e/2 +
f*x/2) - 15*c**3*f) - 282*B*a**3/(15*c**3*f*tan(e/2 + f*x/2)**7 - 75*c**3*f*tan(e/2 + f*x/2)**6 + 165*c**3*f*t
an(e/2 + f*x/2)**5 - 225*c**3*f*tan(e/2 + f*x/2)**4 + 225*c**3*f*tan(e/2 + f*x/2)**3 - 165*c**3*f*tan(e/2 + f*
x/2)**2 + 75*c**3*f*tan(e/2 + f*x/2) - 15*c**3*f), Ne(f, 0)), (x*(A + B*sin(e))*(a*sin(e) + a)**3/(-c*sin(e) +
c)**3, True))

Maxima [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1685 vs. \(2 (150) = 300\).

Time = 0.35 (sec) , antiderivative size = 1685, normalized size of antiderivative = 11.01 \[ \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx=\text {Too large to display} \]

[In]

integrate((a+a*sin(f*x+e))^3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^3,x, algorithm="maxima")

[Out]

-2/15*(3*B*a^3*((105*sin(f*x + e)/(cos(f*x + e) + 1) - 189*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 200*sin(f*x +
e)^3/(cos(f*x + e) + 1)^3 - 160*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 75*sin(f*x + e)^5/(cos(f*x + e) + 1)^5
- 15*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 - 24)/(c^3 - 5*c^3*sin(f*x + e)/(cos(f*x + e) + 1) + 11*c^3*sin(f*x +
e)^2/(cos(f*x + e) + 1)^2 - 15*c^3*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 15*c^3*sin(f*x + e)^4/(cos(f*x + e)
+ 1)^4 - 11*c^3*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 5*c^3*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 - c^3*sin(f*x
+ e)^7/(cos(f*x + e) + 1)^7) + 15*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/c^3) + A*a^3*((95*sin(f*x + e)/(cos(
f*x + e) + 1) - 145*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 75*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 15*sin(f*x
+ e)^4/(cos(f*x + e) + 1)^4 - 22)/(c^3 - 5*c^3*sin(f*x + e)/(cos(f*x + e) + 1) + 10*c^3*sin(f*x + e)^2/(cos(f*
x + e) + 1)^2 - 10*c^3*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 5*c^3*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - c^3*s
in(f*x + e)^5/(cos(f*x + e) + 1)^5) + 15*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/c^3) + 3*B*a^3*((95*sin(f*x +
e)/(cos(f*x + e) + 1) - 145*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 75*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 15
*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - 22)/(c^3 - 5*c^3*sin(f*x + e)/(cos(f*x + e) + 1) + 10*c^3*sin(f*x + e)^
2/(cos(f*x + e) + 1)^2 - 10*c^3*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 5*c^3*sin(f*x + e)^4/(cos(f*x + e) + 1)^
4 - c^3*sin(f*x + e)^5/(cos(f*x + e) + 1)^5) + 15*arctan(sin(f*x + e)/(cos(f*x + e) + 1))/c^3) + A*a^3*(20*sin
(f*x + e)/(cos(f*x + e) + 1) - 40*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 30*sin(f*x + e)^3/(cos(f*x + e) + 1)^3
- 15*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - 7)/(c^3 - 5*c^3*sin(f*x + e)/(cos(f*x + e) + 1) + 10*c^3*sin(f*x +
e)^2/(cos(f*x + e) + 1)^2 - 10*c^3*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 5*c^3*sin(f*x + e)^4/(cos(f*x + e) +
1)^4 - c^3*sin(f*x + e)^5/(cos(f*x + e) + 1)^5) - 9*A*a^3*(5*sin(f*x + e)/(cos(f*x + e) + 1) - 5*sin(f*x + e)
^2/(cos(f*x + e) + 1)^2 + 5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 1)/(c^3 - 5*c^3*sin(f*x + e)/(cos(f*x + e) +
1) + 10*c^3*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 10*c^3*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 5*c^3*sin(f*x
+ e)^4/(cos(f*x + e) + 1)^4 - c^3*sin(f*x + e)^5/(cos(f*x + e) + 1)^5) - 3*B*a^3*(5*sin(f*x + e)/(cos(f*x + e)
+ 1) - 5*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 1)/(c^3 - 5*c^3*sin(f*
x + e)/(cos(f*x + e) + 1) + 10*c^3*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 10*c^3*sin(f*x + e)^3/(cos(f*x + e) +
1)^3 + 5*c^3*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - c^3*sin(f*x + e)^5/(cos(f*x + e) + 1)^5) + 6*A*a^3*(5*sin(
f*x + e)/(cos(f*x + e) + 1) - 10*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 1)/(c^3 - 5*c^3*sin(f*x + e)/(cos(f*x +
e) + 1) + 10*c^3*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 10*c^3*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 5*c^3*sin
(f*x + e)^4/(cos(f*x + e) + 1)^4 - c^3*sin(f*x + e)^5/(cos(f*x + e) + 1)^5) + 6*B*a^3*(5*sin(f*x + e)/(cos(f*x
+ e) + 1) - 10*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 1)/(c^3 - 5*c^3*sin(f*x + e)/(cos(f*x + e) + 1) + 10*c^3
*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 10*c^3*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 5*c^3*sin(f*x + e)^4/(cos(
f*x + e) + 1)^4 - c^3*sin(f*x + e)^5/(cos(f*x + e) + 1)^5))/f

Giac [A] (verification not implemented)

none

Time = 0.44 (sec) , antiderivative size = 215, normalized size of antiderivative = 1.41 \[ \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx=\frac {\frac {30 \, B a^{3}}{{\left (\tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 1\right )} c^{3}} - \frac {15 \, {\left (A a^{3} + 6 \, B a^{3}\right )} {\left (f x + e\right )}}{c^{3}} - \frac {4 \, {\left (15 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{4} + 45 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{4} - 30 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{3} - 210 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{3} + 100 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 420 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} - 50 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 270 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 13 \, A a^{3} + 63 \, B a^{3}\right )}}{c^{3} {\left (\tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 1\right )}^{5}}}{15 \, f} \]

[In]

integrate((a+a*sin(f*x+e))^3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^3,x, algorithm="giac")

[Out]

1/15*(30*B*a^3/((tan(1/2*f*x + 1/2*e)^2 + 1)*c^3) - 15*(A*a^3 + 6*B*a^3)*(f*x + e)/c^3 - 4*(15*A*a^3*tan(1/2*f
*x + 1/2*e)^4 + 45*B*a^3*tan(1/2*f*x + 1/2*e)^4 - 30*A*a^3*tan(1/2*f*x + 1/2*e)^3 - 210*B*a^3*tan(1/2*f*x + 1/
2*e)^3 + 100*A*a^3*tan(1/2*f*x + 1/2*e)^2 + 420*B*a^3*tan(1/2*f*x + 1/2*e)^2 - 50*A*a^3*tan(1/2*f*x + 1/2*e) -
270*B*a^3*tan(1/2*f*x + 1/2*e) + 13*A*a^3 + 63*B*a^3)/(c^3*(tan(1/2*f*x + 1/2*e) - 1)^5))/f

Mupad [B] (verification not implemented)

Time = 14.76 (sec) , antiderivative size = 336, normalized size of antiderivative = 2.20 \[ \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^3} \, dx=\frac {\frac {52\,A\,a^3}{15}-\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )\,\left (\frac {40\,A\,a^3}{3}+82\,B\,a^3\right )+\frac {94\,B\,a^3}{5}+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^6\,\left (4\,A\,a^3+12\,B\,a^3\right )-{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^5\,\left (8\,A\,a^3+58\,B\,a^3\right )-{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^3\,\left (\frac {64\,A\,a^3}{3}+148\,B\,a^3\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^4\,\left (\frac {92\,A\,a^3}{3}+134\,B\,a^3\right )+{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^2\,\left (\frac {452\,A\,a^3}{15}+\frac {744\,B\,a^3}{5}\right )}{f\,\left (-c^3\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^7+5\,c^3\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^6-11\,c^3\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^5+15\,c^3\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^4-15\,c^3\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^3+11\,c^3\,{\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )}^2-5\,c^3\,\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )+c^3\right )}-\frac {2\,a^3\,\mathrm {atan}\left (\frac {2\,a^3\,\mathrm {tan}\left (\frac {e}{2}+\frac {f\,x}{2}\right )\,\left (A+6\,B\right )}{2\,A\,a^3+12\,B\,a^3}\right )\,\left (A+6\,B\right )}{c^3\,f} \]

[In]

int(((A + B*sin(e + f*x))*(a + a*sin(e + f*x))^3)/(c - c*sin(e + f*x))^3,x)

[Out]

((52*A*a^3)/15 - tan(e/2 + (f*x)/2)*((40*A*a^3)/3 + 82*B*a^3) + (94*B*a^3)/5 + tan(e/2 + (f*x)/2)^6*(4*A*a^3 +
12*B*a^3) - tan(e/2 + (f*x)/2)^5*(8*A*a^3 + 58*B*a^3) - tan(e/2 + (f*x)/2)^3*((64*A*a^3)/3 + 148*B*a^3) + tan
(e/2 + (f*x)/2)^4*((92*A*a^3)/3 + 134*B*a^3) + tan(e/2 + (f*x)/2)^2*((452*A*a^3)/15 + (744*B*a^3)/5))/(f*(11*c
^3*tan(e/2 + (f*x)/2)^2 - 15*c^3*tan(e/2 + (f*x)/2)^3 + 15*c^3*tan(e/2 + (f*x)/2)^4 - 11*c^3*tan(e/2 + (f*x)/2
)^5 + 5*c^3*tan(e/2 + (f*x)/2)^6 - c^3*tan(e/2 + (f*x)/2)^7 + c^3 - 5*c^3*tan(e/2 + (f*x)/2))) - (2*a^3*atan((
2*a^3*tan(e/2 + (f*x)/2)*(A + 6*B))/(2*A*a^3 + 12*B*a^3))*(A + 6*B))/(c^3*f)